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Abstract 

The performance of maximal ratio combining (MRC) in Rayleigh channels with co-channel interference (CCI) 
is well-known for receive arrays which are co-located. Recent work in network MIMO, edge-excited cells and 
base station collaboration is increasing interest in macrodiversity systems. Hence, in this paper we consider the 
effect of macrodiversity on MRC performance in Rayleigh fading channels with CCI. We consider the uncoded 
symbol error rate (SER) as our performance measure of interest and investigate how different macrodiversity power 
profiles affect SER performance. This is the first analytical work in this area. We derive approximate and exact 
symbol error rate results for M-QAM/BPSK modulations and use the analysis to provide a simple power metric. 
Numerical results, verified by simulations, are used in conjunction with the analysis to gain insight into the effects 
of the link powers on performance. 

Index Terms 

Macrodiversity, MRC, symbol error rate, Rayleigh fading, Network MIMO, CoMR 

I. Introduction 

Maximum ratio combining (MRC) is a well-known linear combining technique that maximizes the 
signal-to-noise ratio (SNR) in noise limited systems [1J. In the presence of co-channel interferers, MRC is 
sub-optimal compared to minimum mean squared error (MMSE) combining. However, MMSE combining 
requires instantaneous channel knowledge of both the desired source and interfering sources. In contrast, 
MRC only requires a knowledge of the desired source and hence is simpler to implement. For this reason, 
there is still interest in MRC processing in the presence of interference. In [2J, MRC is investigated 
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for large systems where it was shown that in the limit as the number of antennas increases, intercell 
interference effects disappear. In [3], a switched MRC/MMSE receiver is proposed where the simplicity 
of MRC is preferred when the interference levels drop below a threshold. Here, MRC performance in 
the presence of small but non-zero interference is important. The performance of MRC systems with co- 
located antenna arrays is well known for Rayleigh fading channels with multiple co-channel interferers flU, 
flU. Recently, interest in distributed combining has grown due to research in cooperative systems, base- 
station collaboration [6, pp. 69], edge-excited cells 0, (HI and network MIMO IfTOl . In the standards, 
distributed processing is part of coordinated multipoint transmission (CoMP) in LTE Advanced. For these 
macrodiversity systems, every link may have a different average SNR since the sources and the receive 
antennas are all in different locations. This variation in SNR makes performance analysis more complex 
and to the best of our knowledge no analytical results are currently available for such systems. 
Hence, in this paper we analyze the symbol-error-rate (SER) of macrodiversity MRC systems. Note that the 
system is not new. Standard MRC processing is considered and so the general form of the receiver output 
and the initial steps in the performance evaluation are well-known. However the macrodiversity layout 
creates a new channel structure which is far more complex than the microdiversity channel. Hence, the 
MRC output has a completely new statistical distribution and a novel, more advanced analysis is required 
for system performance evaluation. In particular, we consider a distributed antenna array performing MRC 
combining for a single antenna desired source in the presence of an arbitrary number of single antenna 
co-channel interferers. The analysis also covers the case where both the desired and interfering sources 
may have multiple antennas. Since the sources and the receive antennas are not co-located, the channels 
are normally independent and so the focus is on independent Rayleigh fading channels where each link 
has a different SNR. In this paper, we evaluate the SER over Rayleigh fading for fixed values of the long 
term link SNRs. Hence, the SER is computed over fast fading while path loss effects and shadowing 
are held constant. Looking at the joint effects of the slow fading (see, for example, ifTTTl . Ifl2l0 would be 
an interesting topic for future work. In the scenario where some sources have multiple antennas, there 
may be spatial correlation in the channels corresponding to the antennas at that source. However, this is 
beyond the scope of the current work where independent channels are considered. We provide specific 
results for BPSK and QPSK modulations, but the analysis can be applied to Af-QAM and a wide range of 
modulations where the SER can be written in terms of an expected value of the Gaussian Q-function and 
Q 2 -function. The general analytical approach follows the techniques in [13]. The novelty in the analysis 
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is the identification of a representation for the interference and noise term in the combiner output and the 
use of this representation in exact SER calculations. We then use the SER results to analyze the effect of 
macrodiversity on MRC performance. 

The rest of the paper is organized as follows. In Sec. II, we give the system model and in Sec. Ill the 
performance results and SER are derived. Sec. IV gives numerical results where the analysis is verified 
by simulation and conclusions are presented in Sec. V. 

II. System Model 

Consider N single- antenna distributed users communicating with hr distributed transmission points 
(TP) [fT51 each with a single receive antenna over an independent flat fading Rayleigh channel. The 
system diagram is given in FigfJJ The received signal is given by 

r = Hs+n, (1) 

where r = (ri, r 2 , . . . , r nR ) T is the C nRXl receive vector, H = (ha,) is the C nRXN channel matrix, 
8 — (si, s 2) • • • j sn) T is the C Nxl signal vector and n — (ni, n 2 , . . . , n nR ) T is the C nRXl additive- white- 
Gaussian-noise (AWGN) vector at the receive antennas such that n ~ CM (0, a 2 I). The signals are 
normalized to be zero-mean, unit power variables so that E {\si\ 2 } = 1 for i — 1, 2, . . . , N. The channel 
matrix, H, has independent zero-mean, complex Gaussian elements such that E{\H ik \ 2 } = P ik . Hence, 
equation (OQ) can be rewritten as 

r = (^oP^j + fi, (2) 

where P = (Pik), -P° 5 is the element-wise square root of P, the operator, o, represents Hadamard 
multiplication and the elements, H w , ik , of H w satisfy H Wyik ~ CM (0, 1) Vi, k. The matrix, P, is the 
global power matrix for the system and for the k th source, an individual power matrix is also defined by 
Pfe = E ^hkh k | = diag (Pi*, P2k, • • • , Pn R k), for k = 1, 2, . . . , N. In the microdiversity case, P k oc /. 
In macrodiversity scenarios, P k is no longer proportional to the identity and these more general power 
matrices make the analysis more complex. Assume, without loss of generality, that user 1 is the desired 
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Fig. 1. System diagram. To reduce clutter, only paths from source 2 are shown. 

user. For the purpose of decoding user 1, ([I]) can be rewritten as 

r = h i .s i + Hs + n 
= h 1 s 1 + i, 



(3) 
(4) 



where hi is the first column of H, H is all columns of H, excluding the first column, meaning H = 
(hi,Hj, and s = (s 2 , . . . ,sjy) T . The n R x 1 vector, i is the interference and noise vector. With MRC 
processing, the output of the combiner is given by lfl"3l 



8l + 



hi hi 
The interference and noise term in © can be written as 

_ hfi _ h " {Hs + n 



(5) 



hfhi 



h^hi 



(6) 



Following the standard approach lfT3l . we develop a conditional Gaussian representation for Z as follows. 
Since H and n are zero-mean Gaussian and independent of hi and s, it follows that Z is also zero-mean 



Gaussian conditioned on hi and s. The conditional variance of Z is given by 

[/if (Hs+n) (i H H H +n H ) h x 



E{\Z\ 2 \h u i}=E 



h*E \ks~s H H H + a 2 /} h! 



h^ 1 h\ 



hi,i >, (V) 



(8) 



(9) 

Hence, since Z is a conditional Gaussian with variance given by ©, it follows that Z has the exact 
representation 



n 



Z = ~ TWl U> ( 10 ) 



'hi (l£*Ph\8k\ a + <T 2 l)h 
h^hi 

where U ~ CM (0, 1). Using this representation in © gives the combiner output in simplified signal plus 
noise form as 



r = si + ^U, (11) 
where X = hfh lt Y = fcf D (5) and £> (5) = £f =2 Ffc N 2 + ct2/ - 

III. Performance Analysis 

A. A Simple SER Analysis 

With the combiner output given by (fTTI) . SERs for many modulations can be obtained using standard 
methodology [fT3l . As an example, for BPSK, we have the SER 



P„ = Pr(-l+^Re(E7)>0). (12) 



E<Q\\i^\>, (13) 



where C} (x) = -A= J^ 00 e a dt is the Gaussian Q-function defined in IfToll . Defining 7 = X 2 F 1 gives the 



BPSK SER as E {Q (v 7 ^) }• Note that in general 7 is a function of s but this dependence is not shown 
for convenience. For BPSK, each element of s has unit modulus and so there is no dependence on s and 
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the SER in <fT3T ) is valid for any values of s. For many modulations iTTVTl , |fT8l , SERs are constructed from 
similar functions of the form 



where / (7) is the probability density function (pdf) of 7. For some modulations, such as BPSK, the 
SER can be given exactly in terms of Wi(a,b, I), whereas for other modulations it will provide an 
approximation. Using integration by parts on (fT5l) gives 



where F (.) is the cumulative distribution function (cdf) of 7. Hence, SER performance for MRC relies 
on the evaluation of (TT6l) which in turn relies on the cdf of 7. 

In the microdiversity case, all the F; matrices are proportional to the identity and 7 reduces to a simplified 
expression, 7 oc x 2 , where x 2 is a chi-squared random variable [fT9l . In the macrodiversity case, this 
reduction does not occur and 7 is proportional, not to a simple chi-squared random variable, but to 
a ratio of powers of correlated quadratic forms. This is the novel analytical challenge posed by the 
macrodiversity scenario. The derivation of the cdf is based on the joint distribution of X, Y. From [|20l , 
the joint distribution of X and Y becomes 




(14) 



(15) 




(16) 



fx,Y (x, y) = Yl & 



riR n R 




Qjx\ \u[y - f^x) u (x) for f3 ik > 



(17) 



i=l ky^i 



-uljgx-y)u(x) for f3 ik < 0, 



where u (x) is the standard unit step function defined as 




(18) 



and 

P^ R - 2 v ik n «- 3 

£ifc — y~t n R 7 \ ' (19) 

0ik = —, (20) 

Vik 

Vik = PilQk — QiPkl, (21) 

Vik = Pa — Pki, (22) 
Q = D(s)P 1 = diag (Q 1; Q 2 , . . . , Q nR ) . (23) 

Note that each term in the summation of (TTTT) has its own region of validity depending on the algebraic 
sign of 0H-. For example, when fi ik > 0, the region of validity becomes the infinite region bounded below 
by the x = and y = %^x curves. The ft®. = condition has been ignored since the case of distributed 
users with a single antenna always yields /3 ik 0. The cdf of 7 is defined by 

F 7 (r) = Pr (7 < r) = Pr f^- < r^j , (24) 

= Pr (X 2 -rY<0), (25) 
fx.y (x, y) dxdy, (26) 

v 

where the domain of integration is defined by V = {x, y : x > 0, y > and x 2 — ry < 0}. In Appendix 
lAl the integral in (|26|) is computed giving 

nR n R 
i=l k^i 

where F ik (r) = F} k (r) for fi ik > and F ik (r) = Ff k (r) for (3 ik < 0, where (r) and (r) are given 
in ([28]) and ([29]) and /3 ife is given in (|2Q|) . 



^ (r) = ^ (l - e-f ) + l-^f (l _ , (^)) , 
(r) = ^ fl - - |L /^ e S|ferfl f V^&a ifc ) + erfl ( i.fZl^)) , (29) 



0ik V / 2/3jfc V — f3ik V V / V 2 V — Afe 
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ft* 



2 1 5Pj! 



(33) 



2 1 bPi 



Uik 



(34) 



where 

1 — PikQi Qi , ^>ik s 

Wifc = 5 , aik = -5- + 77^-, (30a) 

-Til -Til ^Pifc 

and $ (s) = 4= J* e~ t2 dt in (|35l) and (l36l) is the standard error function ifToTl . Furthermore, Q (a;) = 
0.5 ^1 — $ (^gjY For convenience, we expand (3 ik in (|2Q|) and also give the result here as 

_j i_ 

^-ju=2 {.^ku Pitt) \Su\ 

Note that the case of /3 iA = is not considered as this is the case when Pn = P kl , an event which occurs 
with probability zero when the receive antennas are not co-located. As for the cdf, the SER analysis is 
performed separately according to the algebraic sign of /3 ik . Therefore, substituting (|271 ) into (fT6l) , the 
final result is 

riR n R 

W 1 (a,b,s) = J2J2 P ^ < 32 > 

8=1 ky^i 

where P Sl = P} for ih > and P S l = P 2 for p ik < 0, and P, 1 and P 2 are given in O and (El). 
The results in (|33l) and (1341) are obtained using the following three standard integral identities lfT6ll 



jf (1 - $ (Vqx))& = ^ (^/^~ " - X ) ' for Re (a) > 0; Re (//) < Re (a) , (35) 

/*oo 

/ xe^ x,2 $ (jax) = J ° for Re (//) > 0; Re (» > Re (a 2 ) , (36) 

Jo 2/iy/i — a 2 



00 Ar 

e^'dx = ^_ for g > 0. (37) 
2g 



For multi-level constellations, the values of s affect D (s) and therefore 7. Hence, SER results must 
average (|32b over all possible values of 5. This gives 

W 1 (a,6) = ^W 1 (o,6,i)Pr(i), (38) 
s 

where (138b may be an exact or approximate SER result, the summation is over all possible 5 and Pr (5) 
is the probability of a particular s value. Finally, for BPSK modulation, the SER in (fT3b becomes 

P, = Wi (1,2,5). (39) 

B. Extended SER Analysis 

For M-QAM, first order SER approximations can be found via expressions of the form in (fT4b . Exact 
results involve expectation over the Q 2 (-) function in addition to ([141) . As an example, consider 4-QAM 
where the SER is given by 

P s = Pr f — l - + -^Rc (U) > or — l - + -Lira (C/) > J (40) 
V V2 v 7 ^ V2 v 7 ^ / 

— \ 2 
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P. = l-Pr[Re(tO<y^J (41) 

= i-p{(i-g( v ^)) 2 } (42) 

= 2£{Q( v / 7 -)}-£{Q 2 ( v / 7 -)}. (43) 

Here, the 2E {Q (\/7) } = VVi (2, 1, s) term in (|43l is a good approximation to P s [[T31 and the remaining 
term, P {Q 2 (-v/7) }, makes only a small adjustment. However, in other variations of M-QAM modulation 
schemes the contribution from Q 2 (.) is not negligible lfT3~l . Therefore, for general M-QAM, the exact 
SER is useful and this can be written in terms of Wi (a,b,s) and W2 (a, 6,5) = -E {aQ 2 (-v/67)}. The 
first expectation is found in d32l) . The second expectation can be derived as follows. Let, 



W 2 (a, 6, 5) = a J™ Q 2 [y/b^j f (7) d 7 . (44) 
Using integration by parts on (144b gives 



/2" /"°° 2 f w 2 \ 

W 2 {a,b,s) = a\j- J e~Q (w) F 7 \— j dw. (45) 
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In order to facilitate our analysis we need two fundamental probability integrals. Therefore, we derive both 
integrals in Appendix |B] along with their regions of convergence, since they may have applications in other 
communication problems. Note that similar results may be found in lfl4ll . but these are for restricted ranges 
of the parameter values. The macrodiversity integrals require a wider range of values and the analysis in 
Appendix |B] enables us to evaluate both the integral values and the precise region of validity. As for the 
simple SER analysis, the extended analysis is also performed separately according to the algebraic sign 
of (3ik- Therefore, substituting (1271 ) into (|45l ), the final result is derived in Appendix |B] as 

W 2 (a,b,~s) = J2J2 P ^ (46) 

where P Sik = P^. k for f3 ik > and P Stk = P s 2 . fc for {3 ik < 0, where P}. k and P^ k are given in d47j and (|48l) . 
respectively. Hence, the exact SERs are computable using d32l) and (l46l) for any Af-QAM modulation. 
As in Sec. IIII-Al for multi-level constellations the SER results depend on W\ (a, b) and W2 (a, b) results 
where Wi (a, b) is given in (|38l ) and 

W 2 {a,b) = ^W 2 (a,6,s)Pr(s). (49) 
s 

For QPSK modulation the SER in (@3]> becomes 

P t = W! {2,l,i) -W 2 {1,1,8). (50) 

C. A Simple Power Metric 

The SER and any other performance metrics are functions of the power matrices P 1 ,P 2 , . . . ,Pn- 
Although (|32l) and (l46l) give the exact SER as a function of these powers, the result is complex and 
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does not offer any simple insights into the relationship between performance and the powers. Hence, we 
consider © and © which give the mean SINR of the combiner as 

2 



m P ^E{ — ^ ^ ^— >, (51) 

where mp is a performance metric based on the link powers and we have used E {|sj| 2 } = 1. Exact 
evaluation of (|5D) is possible but it is rather involved and produces complex expressions. Hence, we prefer 
the compact approximation based on the first order delta method, similar to the Laplace approximation 
ED, given by 



hi 





{(Khl) 


■! 


\ 




[hi (j2k=2^khk 




Ifcl} 



mp = — p -. — r r-. (52) 



Using established results for the moments of quadratic forms lfT9l pp. 119], we obtain 

Tr(P 1 ) 2 + Tr(P 2 ) 
m P = — ^— v . (53) 

From fl53l) . we observe that while a Pi matrix with large entries boosts the numerator, hence improving 
performance, it also interacts with the interferers in the denominator. Since MRC is based on weighting 
the strongest signal, the most advantageous interference profile is for the stronger interferers to line up 
with the weaker desired signals and vice-versa. This intuitive result is precisely captured by (1531) which 
increases with Tr (Pi) and also increases as Tr (PiPi) decreases, for i — 2, 3, . . . , N. 
Although mp captures some of the important relationships between performance and the power matrices, 
it is not always an accurate predictor of performance. As the SINR grows, the mean becomes further from 
the lower tail which governs error rate performance. Hence, we expect the mean to carry less information 
about SER at high SINR. This is discussed in more detail in Sec. [IV] 

D. Remarks on Systems with Multiple Co-located Receive Antennas at TPs 

The analysis in Sec. Ill] is restricted to situations where the TPs have a single antenna each. However, if 
the receiver, for example has two co-located antennas at any TR the system analysis still can be handled 
by the same method, but will result in a different joint distribution for (fTTT) . Thus, every new scenario for 
co-located antennas gives a new joint distribution and in turn this gives a different error rate expression. 
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A pragmatic solution is to use a perturbation approach. If P a = P rl (corresponding to receive antennas i 
and r of the desired user, being co-located at i th TP) then we can use P a and P ix + e for the two powers 
where e is a small perturbation. This approach provides stable and accurate results as will be shown in 
Sec. 13 

IV. Numerical and Simulation Results 

For the numerical results, we consider a system with three distributed receive antennas and also a larger 
system with 6 receive antennas deployed in three sets of co-located pairs. Hence, there are three positions 
at which one or two antennas are deployed and these are refereed to as locations. Note that the number of 
interferers in the system is irrelevant, since, from (TTOl ), their effect is governed by J2k=2^k \ s k\ 2 - Hence, 
one interferer with a power matrix equal to ^2 k=2 Pk \ s k\ 2 i s equivalent to N — 1 interferers with power 

l 2 i 2 

matrices P 2 |s 2 | , • • • , Pn \sn\ ■ Hence, we consider a single interferer throughout. In this section we 
consider BPSK and 4-QAM results where |sj| 2 = 1 Vz. Hence, for both systems, we parameterize the 
performance by three parameters which are independent of the transmit symbols. The average received 
signal to noise ratio is defined by p = Tr(Pi) /riRa 2 . The total signal to interference ratio is defined by 
q = Tr(Pi) /Tr(P 2 ). The spread of the signal power across the three locations is assumed to follow an 
exponential profile, as in 11221 . so that a range of possibilities can be covered with only one parameter. 
The exponential profile is defined by 

P lk = K k (a)a l -\ (54) 

for receive location i and source k where 

K k {a) =Tr(P fe )/(l + a + a 2 ), A; = 1,2, (55) 

and a > is the parameter controlling the uniformity of the powers across the antennas. Note that as 
a — >■ the received power is dominant at the first location, as a becomes large (a 1) the third location 
is dominant and as a — > 1 there is an even spread, as in the standard microdiversity scenario. In Figs. [2H3] 
we show SER results for the ten scenarios (S1-S10) given in Table HI Note that an error floor occurs as 
p — > oo (a 2 — > 0) for fixed q. The value of the error floor is obtained by letting a 2 — > in (|32l . In Table 
H] we report the values of nip and the error floor for the ten scenarios considered. Note that m p is given 
for a cr 2 value corresponding to p = 20 dB. 
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TABLE I 

Parameters for Figures[2]and[3] 







Decay Parameter 






Sc. No. 




Desired 


Interfering 


rap 

(dB) 


Err. Floor 


SI 

S2 


1 
1 


65 

a = — 

ftp: 

65 

(T = — 

65 

a = 1 


a = i 

DO 

a = 1 


3.06 
7.68 


1.36e-l 
6.26e-2 


S3 


1 


a = 65 


28.64 


1.80e-3 


S4 


1 


a = 1 


5.97 


2.49e-2 


S5 


1 


a = 1 


a = ±- 

65 


5.97 


2.76e-2 


S6 
S7 


10 
10 


" 65 
" 65 


" 65 

a = 1 


12.93 
17.30 


1.42e-2 
4.90e-3 


S8 


10 


" 65 


a = 65 


27.62 


1.68e-4 


S9 


10 


a = 1 


a = 1 


15.60 


1.21e-4 


S10 


10 


a = 1 


a — 65 


15.60 


2.57e-4 


TABLE B 
Parameters for FigureO 






Decay Parameter 






Sc. No. 


? 


Desired 


Interfering 


mp 

(dB) 


Err. Floor 


Sll 
S12 


30 
30 


" 65 

f> = — 

u 65 

a = — 

65 

a = 1 


" 65 

a = 1 


17.42 
21.34 


1.54e-2 
5.20e-3 


S13 


30 


a = 65 


27.68 


1.99e-4 


S14 


30 


a = 1 


19.65 


7.68e-5 


S15 


30 


a = 1 


^ = 65 


19.64 


1.72e-4 



Figures |2| and [3] verify the analytical results in (|32l) for BPSK modulation with simulations and also 
explore the effect of different power profiles. In Fig.[2l a low SIR is considered with q = 1. Here, SI is the 
worst case since the desired signal profile is aligned with the interferer and the profile is rapidly decaying 
giving little diversity. S3 is the best since the profiles are opposing and the best desired signal aligns with 
the weakest interference. Since Fig. [2] has a low SIR the major impact on performance is caused by the 
presence or absence of a high SIR or low SIR at each antenna. In Fig. [3j the same power profiles are 
considered but at higher SIR, q = 10, the order is changed. S6 is still the worst as this scenario has high 
interference at all antennas and little diversity. In contrast S8 is no longer the best with S9 now giving 
better performance. Note that S9 has greater diversity with an even spread of power across the antennas 
and this becomes more important at high SIR. 

Another comparison between scenarios can be seen in Table HI Note that in Fig. [2] the ordering based on 
mp correctly identifies the best and worst scenarios whereas in Fig. [3] the m p metric suggest that S8 is 
best whereas S9 is better. The m p metric gives some intuition about macrodiversity MRC performance, 
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TABLE III 
Parameters for Figure[5] 







Decay Parameter 






Sc. No. 




Desired 


Interfering 


mp 

(dB) 


Err. Floor 


S16 
S17 
S18 
S19 
S20 


20 
20 
20 
20 
20 


65 

a = — 

65 

(T = — 

" 65 

a = 1 
a = 1 


a = i 

65 

a = 1 
a = 65 
a = 1 

" 65 


17.57 
21.69 
29.32 
20.57 
19.96 


1.50e-3 
1.61e-4 
5.51e-7 
1.54e-7 
1.04e-6 



si 

kt~ — •- . . . . y 




• • • — • ii 

S2 . 

' . . . ^ n 




S5 




S4 . 

*\ S3 


Simulation 




• Analytical (exact) 





p[dB] 

Fig. 2. Analytical and simulated SER values for a MRC receiver with BPSK modulation in flat Rayleigh fading for scenarios S1-S5 with 
parameters: na = 3 and ? = 1. 



especially at low SIR, but it doesn't accurately capture diversity effects (seen in the lower tail of the 
combiner output) which are needed for accurate performance prediction. 

In Fig. the same power profiles are considered for QPSK transmission. Here, the exact results from 
Sec. IIII-BI are verified by simulation. In particular, the SER expression in (1431) for QPSK modulation is 
used along with (T32l) and (l46l) . The relative performance provided by the 5 scenarios is the same as in 
Fig. [3] except that the cross over of S13 and S15 in Fig. |4] (equivalent to the cross over of S8 and S10 in 
Fig. [3]) does not occur until SNR > 30 dB. 

Finally, in Fig. |5] we consider the six antenna receiver where antennas 1,2 are co-located, antennas 3,4 
are co-located elsewhere and antennas 5,6 are also co-located and separated from antennas 1-4. Here, the 
long term receive SNR of a source at antennas 1 and 2 will be the same. Hence, we use the perturbation 
approach of Sec. lIII-Dl to obtain results. Fig. [5] validates the perturbation approach by simulation and shows 
a large performance improvement relative to Fig. |4] due to the increased number of antennas. Again, the 
results due to the five scenarios follow the same order as in Figs. [3] and HI Note that when a = 1 for 



15 



10° 
10" 1 

a 10 

CO 

10- 3 
io- 4 

-5 5 10 15 20 25 30 

p[dB] 

Fig. 3. Analytical and simulated SER values for a MRC receiver with BPSK modulation in flat Rayleigh fading for scenarios S6-S10 with 
parameters: ur = 3 and ? = 10. 
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Fig. 4. Analytical and simulated SER values for a MRC receiver with QPSK modulation in flat Rayleigh fading for scenarios S11-S15 
with parameters: ur = 3 and c = 30. 

both desired and interfering sources, the system layout is microdiversity. Hence, scenarios S4, S9, S14 
and S19 provide microdiversity results. 

V. Conclusion 

Exact SER results are derived for BPSK and M-QAM modulations in a Rayleigh fading macrodiversity 
system employing MRC. The results have applications to several systems of current interest in commu- 
nications including network MIMO and cooperative communications. The analysis is used to study the 
effects of the macrodiversity power profiles on MRC performance. It is shown that simple power metrics 
may capture several features of MRC performance but the impact of diversity in a distributed system is 
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Fig. 5. Analytical and simulated SER values for a MRC receiver with QPSK modulation in flat Rayleigh fading for scenarios S16-S20 
with parameters: hr = 6 and ? = 20. 



important at realistic SINR values. Here, the exact results are necessary to provide an accurate performance 
measure. In general, performance improves as the desired signal dominates the interferer at some antennas 
and as the desired power is spread more evenly over the receive antennas. The exact balance between 
these two key features is difficult to obtain in a simple form but is provided by the exact solutions given. 

Appendix A 
Calculation of the cdf of 7 

Since each term in the summation of (fTTI) depends on the algebraic sign of /3 ik , the final cdf has two 
parts as below 

n R n R 

where F ik (r) = F/ fc (r) for (5 ik > and F ik (r) = Ff k (r) for /3 ik < 0. In subsection IA-AI we derive F^ k (r) 



followed by the derivation of Ff k (r) in subsection IA-BI 



A. Derivation of Fl k (r) 

From the joint pdf in (fTTT ), when /3 ik > 0, F/ fe (r) is given by 



fx,Y (x, y) dxdy, 



(57) 
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where T\ = |x, y : x, y > 0, y — j^x > 0, y — ^ > j. By using standard methods for 2-D integrals we 



arrive at 

oo 



*4 (r) = ~ & T / ' e-^eM^) dydx . (58) 



The final result then becomes 



where 



Flk ,., = ^ (, _ e -t ) + |_ /£ e 5| (l _ . (v ^ aa) ) , 



1 PikQi ,,_r\ \ 

u ik = , (60a) 

-Til 

°" = flT + W <60b) 



The expression in (1281) follows using standard methods of integration in (1581) and employing the 
following integral identity |[T6ll where necessary: 



OO 1 / — 

2 , 1 In 



e-P* dx = -J- (1-$ (y/?a)J . (61) 

B. Derivation of Ff k (r) 

From the joint pdf in (flTT) . when j3 ik < 0, Ff k (r) is given by 

f * (r) = /X /x - y (x> y) (62) 

where J-2 = y '■ x,y > 0,y — j^x < 0, y — ^ > o|. By using standard methods for 2-D integrals we 
arrive at 



Pil / Pil 



Fl (r) = J ' J 2 " e-^e-^ y -^)dydx. (63) 
The final result then becomes 

Fl (r) = ^ f 1 - e"f 1 - & /^e^lferfi f v 73 ^^) + erfi f ^/4^Y| , (64) 
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where 



erfi (x) 



$ (jx) 



(65) 



J 



The function erfi (.) is the error function with a complex argument defined in IfTBll . Note that the square 
roots appearing in (l29l are the positive square root of /3^. The expression in (|29l) follows using standard 
methods of integration and employing the following integral identity IfTBI where necessary: 



The integral in (1451) is required for the exact SER analysis. Substituting F y [w 2 /b) from (|T7l) into (|45l> 
gives two new integrals involving Fl k (w 2 /b) or F 2 k (w 2 /b), which are given in (|28l) and (|29|) . These two 
integrals can be written in terms of known functions and two fundamental probability integrals that we 
denote I\ (a,/3) and I 2 (ot,/3). These integrals are computed below. 

A. Integral Form I 
Consider the integral, 




(66) 



Appendix B 



Derivation of the exact SER 




(67) 



Applying the integral forms of Q(.) and $(.) gives 




(68) 



Using the substitutions, t\ = rcosO and t 2 = rsin^, the integral then becomes 




xre 




(69) 
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where tan0 = a\/~2, r\ = x/y/2cos9 and r 2 = ax/ sin 9. Using standard methods of integration with 
some simplifications we obtain 

• raV2 J f - 2 /-on/2 j f 

A (a 8) = / - — / - (70) 

U ' 2tt8J t 2 - 1-28 2tt/3 7 [a 2 - B)t 2 - a 2 

Defining 

f a ^ dt 

IM = L W^Tp- (71) 

f 01 ^ 2 dt 

h2(a,B)= j-z flW9 - , (72) 



o (a 2 - 3) t 2 - a 2 ' 



allows d70l) to be rewritten as 



2 

h (a ' /?) = 2^ hl (a ' /?) " 2^ Jl2 ( "' /?) ' (73) 



The integral in (1711) and (1721) can be solved in closed form to give 

1 , , _j / «V2 



I u (a, /3) = tanh" 1 - (74) 



and 



7 12 {a, 8)={ " (75) 
^=tan~ 1 (j2 (6 -a 2 )) otherwise. 

a^J fi-a 2 \ v / 

Note that some intermediate steps in the derivation show that 1 + 28 > 2a 2 is required for the existence 
of (167T ). This constraint is satisfied by the current problem. This can easily be seen by substituting the 
arguments of both I x (., .) functions in (1481) in to 1 + 2/3 > 2a 2 followed by simplifications using (1301) . 

B. Integral Form II 
Consider the integral, 

/■oo 

l 2 (a,0)= xe^Q(x)(l-^(ax))dx. (76) 
Jo 

Applying the integral forms of Q{.) and $(.) we obtain 

r\ POO POO POO 

I 2 (a,8) = - / / xe^-^dhdhdx. (77) 
vr Jo J * J ax 

V2 



20 



Following the same procedure as in Appendix IB-AI and with some simplifications we arrive at 



f ,. oP_ r+ dd 1 r' 2 do i_ 

2{a ^ } ~ 2ttP J a 2 -/3sm 2 6 + 2tt/3 l-2/3cos 2 c? 4/3' <;) 
where tan</> = y2a. Making another substitution as t — tan# in (1781) gives 

, _ a 2 r^ a dt i r 00 rft i_ 9 



Defining 



M«,fl=/ 7-5 ^ ■ _„ (80) 



o (a 2 - /3) t 2 + a 

^22 («,/?)=/ , (81) 



allows (|79|) to be rewritten as 



h (a ' /?) = 2^/ 21 ( "' /?) + 2^ /22 ( "' /?) " 4^ (82) 



The integrals in (1801) and (1811) can be solved in closed form to give 



V2 ^2 



a 2 = /3 

/21 (a, /?)=<( " (83) 
/ tan" 1 (J2(a 2 -/3) ) otherwise 



and 



i 



1 = 2/? 

/, 2 (n. /) = < VZQ (84) 

7ltT coth_1 (^r) otherwise. 
As for ii (a, /3) there is an associated region of validity, 1 + 2a 2 > 2/3, which is satisfied by the problem. 
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